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soiEivjor 

The unstable 'burning of a bipropellant rocket combustion 
chamber Is investigated end a study made of the requirements for an 
automatic closed loop control circuit to stabilize the motor. 

Coe bipropellent combustion chamber equations developed ty 
Dr. L. Croceo^ are utilised as tho analytical description of the 

rocket motor burning phenomena. Squat lone similar to those developed 

(o) 

by Dr. H. S. Talon aro used for the oxidises* and fuel supply systems 
and tho two dosod loop stabilising circuits. 

Ihe stability or instability of the system is demonstrated by 
the use of a special plotting diagram in the complex plane suggested 
by it. Satche as a moans of handling ay 3 tome with time lag, and 
developed for thl 3 use 'Ey H. S. Taien. This involves separating the 
transfer function into two parts. In tho complex plane the first 
portion of tha transfer function, the exponential variable containing 
the time lag, plots as a unit circle as the complex variable p is 
made to take a contour enclosing the positive half of tho p-plane. 

If tho remaining portion of the transfer function intersects this 
unit circle, the rocket motor can be unstable for large reduced time 
lag; if it does not intersect the unit circle, the system is 
generally stable, although the roots of the exponential coefficient 
in the positive half of the complex plane must be investigated. -This 
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It fct rocai re: sent cm bo conveniently ceon dished by the -‘-.id of 

a ysuist diayran. 

The o attritions for the feedback circuit arc devo3.o_.-ad end the 
oxidisor end fuel transfer function ren'irenonts are dote rained. 

Tvo cases of stable conbucticn cad f/o cases of unstable 
combustion m*e analyzed. One unstable ease is stabilized by the 
ad. it ion of a feedback circuit. 
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SYMBOLS AJTD JKETTITIOnS 



tine 

instantaneous value of the time lag 

pressure exponent of pressure dependence of the 
processes taking place during the tino lag 

instantaneous pressure in the combustion chamber 

X i ^ n root of the charActeristic equation 
vith the reduced tine as the independent variable 

pressure in the combustion chamber in steady 
operation 

(p - p)/p ~ fractional .variation of pressure 
in the combustion chamber 

instantaneous rate of infection, burning, and 
ejection of propellants 

common value of the above rates of Injection, 
burning and ejection in steady operation 

(*. - 5)/5 *= fractional variation of injection 
rate 

instantaneous mass of gases in the combustion 
chamber 

I'lgf n =» gas residence tine in steady operation 

O „ + f - total residence time of propellants 
in steady operation 

absolute temperature of combustion gases 



vi 

reduced amplification coefficient 
reduced angular fee prency 

instantaneous nixturo ratio 
nlxture ratio in stably operation 

coefficient related to the variation of 

temperature with fixture ratio 

non-dimensional reservoir capacity for the 

oxidizer line 

non-dinonsional reservoir edacity for the 

fuel line 

(r - 1)/ 2 (r+ 1) i- coefficient rs. mating 
the deviation from unity of the nixturo ratio 

subscript indicating injector end 

subscript indicating ejector end 
subscript indicating oxidiser 
subscript indicating fuel 
t/0 » reduced tin© 

^Iq = reduced tine lag 

S 

constant relating the mss flow rate to the 
pressure in the feed system lines 

^A/SApA^ 0^ = inertia parameter in the line 

p/2 A p as pressure drop parameter 

ratio of functions in (d /ds) describing tlie 
response of the feedback circuits 

H + 2K - Ke” p = 



a convenient grouping of toms 
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A* , B * , D», L», M* 


= conples constants used in the trial 
solution 


a(p) 


s= system transfer function 


s x (p)t g 2 (p) 


= components of O(p) used in the Tslen-lntche 
35iagraia 


Ag 0® o 


« convenient system constants defined in 

equations (16) through (21) 


L, M 


« convenient grouping of system const nt s >nd 
variables defined to equations (22) and (3) 


$ 


- F 0 (p)L - Ff(p)H = A»i/p * Sq * + jjj 

a convcsnient symbol for the feedback circuit 
transfer function 



i. i::.:r.::ucTion 



Jacket motors often exhibit rough con oust ion cither as 
chugging (low frequency oscillation "below 100 cycles per second) 
or screening (high frequency oscillations several hundreds of cycles 
per second). 'The first of these phenomena has been investigated 
sufficiently that several authors have developed equations describing 
analytically the oscillation processes in rocket motors^. Those 
aquations which are utilized in this paper are deduced by L. Grocco 
under the following assumptions: 

1. The processes of transfowilng the propellents into high 
temperature gases require a certain tine called the tine leg, / f ; 
it is the tine after which the liquids nay he considered suddenly 
transformed into liot gases. She rate of nixing of the propellants 
depends essentially on the injection system used and recirculation 
effects, but probably not on chamber pressure, lie conditions, 
droplet size, temperature, vaporization and eventual chemical trans- 
formation of the propellents after mixing depend considerably on 
chamber pressure. 

2. The time lag is considered the same for all particles of the 
propellants; the residence tine is considered the sene for all 
particles, l.e., all particles travel from the injoctcr end where 
burning occurs to the exhaust end in the same residence time interval. 






% 



5. 



Ci.-. 1r 



to be: 




f (p) dt ' 



constant 



( 1 ) 



whore p io the chamber pressure end f(p) Is tho tine rote of 
preparation for combustion. 



4. Pressure variations arc sufficiently slow throughout the 

co .'auction chamber that a. unifora jjrossuro can be used, p — p. ( t / . 

5. 'Xerrparatura of the upases is assumed to vary with tine and space 
but not with pressure oscillations. 

6. Analysis of the stabilising circuit itself and its oscill tioas 
affecting tho racket notor are not investigated. 

L, Orocco^^ lias demonstrated that a rocket combustion chnnbor 
is intrinsically unstable, i.e., the burning phenomena itself can 
generate oscillations o von without variation in the injection rate, 
-his occurs because any oscillations of the pressure will cause an 
oscillation in the tine leg by equation (1). ?2hua if the time lag is 
oscillating, it is immediately seen that when the tine lag is 
decreasing, the burning of the particles that were injected later will 
catch up dth the burning of those injected earlier; this increases 
the rate of burning with ros-oct to tho average, diailnrly, the 
opposite happens for an increasing tine lag. 2hu3, if the increased 
or decreased rate of burning should coincide with a pressure 
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incroase or doorcase, 99lf-excito& oscillations can exist resulting 
in unstable burning. 

For the bipropellant case assumption 2 that O is the 
residence tine of all particles means that each particle bums at 
the injector end and travels to the exhaust end car lying with it the 
temperature developed at the combustion instant. Shna, the 
temperature at the exhaust, 2 ?^q, at the tins t corresponds to the 
mixture ratio of the propellants injected at the earlier tine 
t-0^=t-'T' - 0 , Shis can be put in terns of the reduced 
tine by dividing by Og giving: * - - 1. Soorefore, in effect, 

we have two reduced tine lags, & and 1.0. It is the second tine 
lag that introduces additional complications in the formulas for 
the bipropellant case by appearing as an exponential tom in g 2 (p). 

She problem of stabilising a combustion chamber is simplified 
by the use of the Tsiea-Satche Diagram. If the system transfer 
function, G(p), is separated into g^(p) = e” l ’ and gg(p) = all 
remaining terns in G(p), then as the complex variable, p « A * i<*>, 
is made to take the contour enclosing the positive half of the 
p-plane, stability can be conveniently determined if the gg(p) curve 
lies completely dear of the unit circle, g^(p). A llyqp list Diagram 
of the denominator of g 2 (p) also is required to determine if there 
are any poles in the positive half plane. 

For an unstable case in which gg(p) does intersect the unit 
circle, g^(p)» stability can be insured by designing a feedback 



circuit transfer function which will cause the go(p) curve to 
so ve clear of the g^(p) curve, and still meet the requirements 
of the llyquist Diagram concerning poles in the positive half TJlaae. 

For a bipropollaat rocket the temperature of the combustion 
gases is a function of the local mixture ratio. As the mixture 
ratio operates away from the steady-state value, r, the fractional 
variation of the temperature varies as: 




or, is- -L ^ 

K ~ * T % dr 



Shis defines K which turns out to be an interesting parameter 
of the system. K nay be either positive or negative, although, 
as pointed out by L. Crocco, K is generally positive. 



II. TCSr-JLOIJrrT OF 3.JJATI0H3 

Under the assumptions given above, the equation of the combustion 
chamber for a bipropellant rocket was developed by L. Groccos (3eo 
Figure 1 for schematic diagram of system). 

3z + Z ~ J V/ SH -1*0-^%] + ?-*<*•- = 

+ ~h h (?- <f ! s ) (2) 

Following H. 3. Tsien and considering the propellant to be fed 
by a centrifugal pump run at constant speed the following equations 
for the bipropollant feed system are developed* 

-i(p t> +iKj^+ = m> + Jo £(«.+ ( 3 ) 

(4 > 

Bquations (2), (3) and (4) specify the transient behavior of 
the fuel and oxidizer supply systems coupled with the combustion 
chamber when the reservoir capacity in the feed line, K 0 and Kg, 
are specified. The closed loop stabilization circuit can bo 
completed by requiring Kq and Kf to be a function of 4* . Tills 



means that the reservoir capacity responds to pressures from the 



combustion chnriber through a suitable amplifier and serve system. 
2Sma Kq end Kf are assumed to satisfy a pair of int ogre-di f f e yon bird 



equations. 



F„ 




- Ko 



(5) 



If <^>4’ “*f 



( 6 ) 



where and ? 0 correspond to a ratio of functions (generally 
polynomials) in d/ds. 

2rying o. solution of the fora: 



0 = A'^ 1 




= «'e^ z 




= 2>'e^ 


K. = 


/ < 1 yZ 
L e r 


Kf 


. /V'e'* 


(7) 


where now p = A 


4- i tO 


is complex, as 


are A*, B*, 


D», L* and h«. 



Sliainaticg the comon factor e 3 , equations (2), (3), (4), (5) 
and (6) represent a solution of the system if the five homogeneous 
first degree equations below are sirnltaaeously satisfied: 



(1 - n + p + ne” ^ )A*-(g + )r)e ^ B* + (g • po ‘ D* *= o 



P 0 A* + [l + J o p ] B * * [<^ p 0 ^i)p + J oP J ]a>* « o 



(9) 



PfA* -t* £l +^(Pj+i') * JfP J D * + |*(Pjf*a')p + JfP 2 * 0 



( 10 ) 



( 11 ) 
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? 0 (p) A* - = 0 



S'f (p) a* “ n* = o 



( 12 ) 



Since they are homogeneous equations, they can be sinul taneoualy 
satisfied if, and only if, the following determinant of the 
coefficients of the above equations vanishes: 



■p + (>- >V + „ 

7>f Sp P Q-i)e A o o 

Po [l+z(Po+i)+J 0 pi o [£ Gi+ i)f+J o p] o 

h(f>) o o -/ C 

Ff(i>) ° o o - / 



where g •= H + 2K » Se”^ has been used for convenience. 

Expanding the determinant gives equation (14) below; if p 
satisfies this equation, our assumed solution, equation (7), will 
be a possible solution of the system. 
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+0-H)J[l +■ K ( P f * y +- Tf?][ l + -< ^oH) + Jo f 3 + 



c '^W/+*(P 4 + i)][/+i(P° + i>] + Po( 9 + < p f + i)]+ 
-P,(yi>[i+i( p -+*>]! + feT*l»%[i+i(ro + i)] + 

" -^Pf^'i) +?f£(9 + £>j + ^ e f ^^ T f + 

u-uim+uUn-* 1 ***'!} ^ 



■^ 5 e ^ f~f(p) ( , j--t'> ( ->f To) +• 

Fi( 1 ,)f[%+ i Ll[«(P°+k)][/+£( p f + ± ) J} f 



?€ s * F(f) (yi) {L±(Z+m I + 



f(t Sp Fo(p (<f^t)( To j+j - G(fi 



(/V 



O 
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for the system to bo ataxic, the roots p of equation (14) 

should not have real parts. The problem of stabilizing the rocket 

motor thus becomes one of desiring the closed loop transfer 

function, T 0 (p) •> F^(p), so that this criterion is satisfied. This 

condition can bo investigated most easily by use of the Tsien-Satche 

'.Diagram and a Hyquist Diagram. 

If the complex variable p is made to take a contour a3 in 

figure 3 enclosing the positive half of the p-plsne from -• O to ■{* cO 

and around a circle in the positive direction with infinite radius, 

each root of CJ(p) = 0 (aquation 14) in the positive half jfLane will 

cause O(p) to make a complete revolution about the origin. Separating 

(J(p) into two parts, g^(p) and g ? (p), by dividing by the coefficient 

_<f p 

of the exponential tom, e , gives* 



= },(f> - Up) 



_ c -^_ f 

(V+Strhp* c+cet*] +(i(p>Uff(p )M \\ ° 5> 



where 



S 1 + £ (Pj. <• |s) 



A 



(16) 
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b s i + i < p 0 + t) 



( 1 ?) 



G = n +(P 0 +P f ) (£t* ^ + £ + £0 * (P 0 -Pf ) <H + 2K + £ i- £ ) + 



Vf (« "*• ^2. ) ^ * 4 *. ^ 






( 18 ) 



D = [(Jf+J 0 ) (tti-^) + (JfP 0 *!*J 0 P f ) (| + £ ) + (Jj^WoV < H + SR?/ & 9 > 



H s (P* - P 0 ) (K + £ ) 



( 30 ) 



Sr -K (J f P 0 - J 0 P f ) 



( 31 ) 



a = 



£(H + 2K - &) [(A - 1) (Bp + JqP 3 ) + J f (Bp 2 ) + 

Vo p3 3 “ Ke ~ P [< A - 1) (*P + J 0 P 2 ) + Jf(Bp 2 ) * JfJoP 3 ]} 

( 22 ) 



is £ (H + 2K + |) [ (B - 1) (A P + J fX P) + J 0 (^p 2 ) + JfJoP 3 ]* 

- Ke“ p [ (B - 1 ) (Ap + J f p 2 ) 4 J 0 (Ap 3 ) + JfJ 0 p 3 ] J 



( 33 ) 
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ITote that In these definitions. A, P, C, I), R, S, are constants 
of the system, hut that L and II include the variable (p). 

With this arrangement of terns, 0(p) becomes a vector with 
vortex on g-^(p) and the starting point on ggCp). As p makes the 
stated contour of Figure 2, g-^(p) becomes a unit circle and g^(p) 
prescribes various contours dependent on its fora. 

«. £ 

It is also necessary to investigate the coefficient of e 
(this is the denominator of g 0 (p) ao written in equation (15) to 
ascertain if it has any poles or aeros in the positive half plane. 

If the coefficient of has s poles and r seros in the 

right half p-plano, then as p traces a clockwise contour as in 
Figure 2, the coefficient will trace r-s clockwise revolutions. 
Therefore, for stability Spiv) must make r counterclockwise revolu- 
tions around the unit circle. Tims even though the primary condition 
of stability is that gg(p) lie completely clear of the unit circle, 
g^(p), it must also satisfy the conditions dictated by the number 
of poles in the positive half plane of the coefficient of the 
exponential term, e” ^ 
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in. STABILITY OF BI-UOi-jaLAIIS STSfcTlS 
•JITHOUT nSDBACK 



Ito investigate first the coupled feed system and combustion 



chamber without the feedback circuit, set K 0 (x>) sad Kj ( P / equal to 
2ero. This reduces equation (15) to; 



r ,,. f H> + '-«]£A+lft>3[B + T.f>3 1 

l-[r’JoJff-+(D<S<i-f’) 1 k+ ) 



(24) 



where the constant 3 A, 2, 0, 3, 3 and B are defined in equations (16) 
through (21). 

0&3B V (Stable Oase, K ~ 0.2) 

let a = 0.1 j ^ = 1.0; r *= 2.75; S - 0.233; = 

dT^/dr s* 650; K = 0.2; J 0 = 4.0; = 0.6; P Q = 

fc> 

P f = 1.0. 

Substituting these values into equations (16) through (21) gives; 

A = 2.5; B = 2.0; C « 1.65; D = 1.08; B = 0.2; 

3 *s 0.72. 

Thus, equation (24) becomes: 



4500°K; 

0.5; 






- J 



(^o. 9 Xi. 5 '+ o.^J(Zo + fo^ 

[0.32 />N- O of + one:P)f TcTisT^z^)] I 



(36) 



Vow let p = 163 , where ^3 is the reduced angular frequsn cy. 
Iquation (25) re due os to: 



Substituting value3 for <-2> gives the Tsion-Satche Diagram for 
g>>(p) plotted In Figure 3. Checking the denominator of go(p) for 
polos gives the Hyquist Diagram plotted In Figure 4. Tho numerator 
is plotted in Figure 5. 

From the figures it Is evident that this case Is stable since 
the go(p) curve is completely clear of the unit circle and makes no 
loops about it. Also, the denominator in the Hyquist Diagram is 
proven not to have any zeros in the positive half p-plene, since it 
does not make any revolutions about the origin. 

2 -.3 ■ II: (Stable Case, K «= 0.3). 

Let n =• 0.1; = 1.0; ? = 2.75; H = 0.233; T - 4500 °K; 

b 

d? g /dr = 1000; K = .3; J 0 = 4.0; J f = 0.8; = 0.5; 

Pj = 1.0, 

Substituting these values into equations (16) throu^i (21) gives: 

A ss 2.5; B a 2.0; 0 = 1.50; D = 0.36; B » 0.3; 3 = 1.08. 

Letting p = 1^0 and substituting these values into equation (24) 




— 3.2. /S’- 4 ^4-cS)c — o^ 1 ' + S' 



(26) 



- 0 . 32 c 2 V(/of ^Orj z ^)i^ JrQ-LS 



gives: 
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Gfr) = e 






_ ) 



(-3.ZU + iS-44<*)L — /^.^cO 2 --/- 4-.<> ) , 

_ f (27) 

-[o.3zoi z +-^.36 + /-oSe- iu, JcoJ -z (A5'+o-3<t <UJ J][ 



Again, substituting values for ^0 gives the Tsisn-Satcho 
Diagran for g^(p) plotted in Figure 6; the Kyquist Diagran of the 
denominator is plotted in Figure 7. She numerator is the 3'«3 as in 
Casa I (Figure 5). 

Changing K in this second case has resulted in a radical change 
in the Toien-Satche Diagran. The curve g.j(p) now makes two counter- 
clockwise loops about the unit circle, although it do 33 not intersect 
the curve g^(p). nevertheless, this is still a stable case since 
the denominator of g>(p) nakes two clodc-rf.se revolutions about the 
origin in the ITyquist Diagram. 

GA35 III: (Unstable Case, K = 0.2) 

■uet n ^ 0.5 2 cK_ 1 .0; r “ 2.75; H -- 0. *33; 
dTg/dr as 650; S = 0.2; J Q e 4.0; J f = 0.3; P Q = 

^ .£ — 1 .0 • 

Substituting these values into equations (16) through (21) gives: 

A - 2.5; B = 2.0; 0 = 3.65; V = 5.72; B = 0.2; 

3 - 0.73. 

Letting p = izO and substituting those values into equation (34) 



4500°K; 

0 . 5 ; 



gives: 



G(p)^ e 




-[-/i o J x -f- (5^72. -f'O.yze t0 ) <.\e (3 6 T +o.ze~ ctM J] 



-h 10 ? oz )i ~ ^3 • 2 to 1 -f- 2-5 



(98) 



Sub3ti tuning values for and plotting gives the "sion-datche 
Diagram of figure 8 for go(p). She behavior of the denominator of 
g^(p) is shown in the plot of Figure 9. The numerator is plotted in 
Figure 10. 

Inspooting the figures it is seen that the g 0 (p) curve now 
intersects the unit circle and thus the rocket motor can be unstable 
for largo reduced tine lags. This shift to instability is primarily 
the result of the change in the value of a fron 0.1 to 0.5. l Fhns, 
n a f, t the pressure exponent of pressure dependence of the processes 
taking place during the tine lag, T' , i 3 shown to be an important 
system parameter. 

CA52 17 l (Unstable Case, X = 0.3) 

Let n = 0.5; = 1.0; r «* 2.75; H *= 0.233; 5? g = 4500° 

d? g /dr = 1000; X « 0.3; J 0 = 4.0; J f = 0.8; F o = 0.5; 



if e 1.0. 

Substituting these values into equations (16) through (21) gives: 



2.5; B - 2.0; C = 3.50; 35 » 4.50; R = 0.3; 



S * 1.08. 
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Letting p = i «■*> and substituting these values into 
equation (24) gives: 



< 3 (p-e 



_ v'cScJ 



r (-1 zcj 3 + /o- ?tu)c - /3-loi 7 - -h Z-& 



^ < 



/.(o tS^-h (4-.S + /.o<2 4 (-3-5 + 0. je' t ^7' 



(29) 



Substituting values for ^ aad plotting the values of g,,(p) 
gives the ? 3 ien-Satcha Diagram of Figure 11* Fho Kyquist Diagram of 
the denominator is plotted in Figure 12 and the numerator, which is 
the sane os in Caso III, is plotted in Figure 10. 

In this ease changing K from 0.2 to 0.3 did not have the 
spectacular results as in the stable case, although the go(p) curve 
does have a more complicated shape. 
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IY. STABILIiATIOII 0? A BirROrELL.MTT SI STS l 
BT A FHKDHA.CK CIRCUIT 

Turning to the problem of stabilising an unstable burning 
rocket no tor. It 13 obvious that a feedback transfer function, 

I? 0 (p)L -t If (p)M, is required that will move the gg(p) curve 
completely clear of the unit circle for snail values of *0 plus 
satisfying the requirements of the Hyquist Diagram for all values of 
^ • Since the gg(p) curve can only intersect the unit circle in 

-i — i^P 

the region of smll values of a -, the © terns in tho denomi- 
nator of (?) (equation 15) — can be expanded in a power series 

of G<*>) and hitler order terns neglected. Also, define a feed 
system transfer function, f Q (p)L + Ff(p)ll, as a series in powers of x^s 

^ 0 (p)b + If (p)M - $ (p) = a-x/p + 3 d + «lP * (30) 

Substituting these two changes into 0(p), equation (15), and 
letting p = i<-3 gives: 



6(p> 



e 



£ cuj *f (?-* y* }JLA+ Jj- fc^j f & + J+ (iu >? ^ 

)]«- -hC+R(i- 



(31) 



How dear tho denominator of go(p) of teaginary terms by 
multiplying both numerator and denominator by the conjugate of the 
denominator giving: 



— « 




[ )][ h Jo r + ^*3 " + <*.3 £m(/~u; -t-Jj &(/~n) + d&J]( ( '<> j ) 1 ' 

( Jz[nJo^ + f-5+'«? J .][Ct^4 <toJ-[V-frS1-4,-]*-}0~)'i [C+X+ a 0 J*- 



{[X>A 4 J f B + J f J« Q.y, )J [C e.i $ -f Lc+rt+4oj[8A(/~>7)3 

{z[nToJf -5 + [Ci- Z^-t-Oo] - [p -R *■ St a,y j (cuj) x -f [C ^ & 0 1 Z 

• J{0(/-n)-h8A]£c+X+<<oJ -[p-R+5*-<x,lC **> ) (3z) 

C + lc+R+<^Y j 

Shis is the approximated equation for the gg(p) term of the 

system transfer function and holds for snail 60 only. Toms of 
, .3 

order (1 <*> ) and higher order were neglected in this development. 

Bow it i3 necessary to specify the shape of the stabilized curve, 
gg(p), near the unit circle and derive the values of a^, a^, aj, ag 
that will meet these requirements. Let this stabilised curve be: 

gg(p) = W + X(i^ ) + I(i^ ) 2 + 3(1 <** ) 3 (33) 

Equate this to gp(p) in equation (31) and clear by multiplying 
the denominator by equation (33), transposing the numerator terns 
to the left hand side and collecting even and odd powers of (icJ ) 
together: 
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+■ -h Y[c**] + Wa z + )(a,-b Ya 0 + Z«-,}(tu>;V 



[iA + Jj: B + (t-n j] (t'oj^-f £y/[C 4 PT ]4 A 6 (/-h) +vja 0 -f X d-r } = 0 ( 34 ) 



{w[l"§-]+x[)o];T f 4-|-5j-l' > r/;P-^53-f-zr^^Ji-X'a 2 +Ya ( + Za 0 }M 3 + 

[l ft ^.](c«-) ? + $ W[p-K+5j-h?< [(+*] + [T 0 A(i-»)+J{ 8 (/->i)+ 5 a] + 

(35) 

Wa, + X* 0 + Y<i-, J(coo) + Vs/<c_i - 0 

Hots that la this development terns of order (ico )^‘ end higher 
were neglected. The coefficient of each (iw) torn met be independ- 
ently aero for these equations to be satisfied. Therefore, sotting 

1 Q 

the coefficient of each tern, i.e., (i 00 >, (ioJ) . (ioJ) , 

P 

(i <*■> )~’ to zero and solving for the unknowns, a_p a^, a^, gives: 

^ a -l = 0 

Wao + X*. 1 = - AB(l-n) - \1 [c + s’] 

1% + + 7a_i = - [j 0 A(l-n) + JfB (l-n>fr3A] -if [ D-n + gj -x[o -4 pj 

Uag <• 3&i *5* Tsq 4* Sa»>i = - [j 0 A + (1-n)] -w[aJ 0 J^ 4* z' +■ 

- X [d - P. * s] • I [c + 2] 



( 33 ) 



These four equations cea be solved siEultsaeously to give 
explicit relations for a_x» Cg, aq» a, in terns of system constants, 
cad thus S’ (p)X ‘‘r ?£(p)ll, the feedback system transfer function, is 
do to mined. 

Solving for the unknown a^, a^, gives: 

a -l - o 

a 0 * - [^ 76 -^)] - (C + B) 

®1 « -[j 0 A(l-a) + J^B(l-n) + B^vS/ - (D-StS) + (C+H) ( X ) * 
[/u3(l-n)^,] 

®2 c * £ J Q A(l-n) * *4* -wj-f nJ 0 j£ ^ “ 8^ *4* 

[^)][£-£] -[o*h][#-£] 

(37) 

For snail 0 the ay.proxina.ted gp(p) equation (32) can "be used 
with the feedback transfer function equation (30) to plot the 
Tsien-Satche Diagram. For largo ^ the general equation (15) must 
be used in plotting this diagram; also tho Ilyquist Diagram of tile 
denominator of g^(p) oust be investigated. 

To illustrate tho procedure outlined above, consider the unstable 
caso discussed previously where K *5 0.2. Let the required stabilised 
curve of go(p) described as follows: 



(38) 



Sg(p) * W+ X(iw) Y(l^ ) 2 * Z(ia>) 3 « - 2 + 8^ 3 + Bid . 

5hia arbitrary selection of tho g 0 (p) carve Is mdo to place 
the stabilised curve in the approximto position of the known stable 
curves of Cases I sad II. Ihus, H = - 2; X = - 8; I = - 8; C - C. 
Also, A = 2.5; B = 2.0; 0 = 3.65; D« 5.7:3; R= 0.2; 3= 0.72 (other 
values 0.9 Case III). 

Substituting these values into equation (37) gives: a_^ t= 0; 

« - 2.6{ = - 5.84; «= - 1.06. With these values, 

± a - 2.6 + 1.06o0 2 - 5.84 ioO. 

using e<jxation (33) and the above values the following expression 
is obtained for £* (p): 




+ 3+3 - IS31 o ) 2 " -h 
-f/.S3 - 3-25* oe 1 



(39) 



Solving this equation for various values of givas the below 
table of results: 
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Vt\ 


W+ X6v))4Y(<c4V 
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-f.se - t 2-ZS 


-f S ~ Z- c»o 


X 
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-/.xi -i 3.5S 


— /. rz- cz. 6t> 
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+ 0-97 - c £7 


G - z 4-. ° 
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Hhis table shows that as becomes greater than 1/3, the 
values of S?(p) have too great an error. However, fron Figure 13 
it Is apparent that the curve has boon stabilised satisfactorily 
for small . 

It Is now necessary to investigate the behavior of g.^(p) 
using the exact equation (15) which includes the feedback circuit; 
substituting the values determined above into equation (15) results 
in: 



Substituting values for and plotting 6g(p) gives the 
‘JJsien-Satclie Diagram of Figure 14. “The Hyquist Diagram of the 
denoainator is plotted in Figure 15; the numerator is the sane as 
Case III and is plotted in Figure 10. 

Fron figure 13 it is evident that the stabilised curve is 
completely clear of the unit circle for small °° . However, for 
large the addition of the feedback circuit has caused the go(p) 
curve to change in the Tsien-Satche Diagram fron no loops about the 
unit circle in Case III to one ’./ith two counterclockwise revolutions 
about g-^(p) in the stabilised curve of Figure 14. nevertheless. 




Sp \ - /3.2oa 2 -f 2.6" -f c (~ ?■ lol 3 -b /o.f OJ) 






( 40 ) 
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the system is still st dole since the Eyquist Piagran shows that 
the denominator of g„(p) makes two clockwise rovolution9 about the 
origin. 

Thus an unstable burning rocket motor has been stabilised by 
an appropriate closed loop feedback system. The circuit chosen 
to stabilise Case III was a combination of three terns of a 
differentiating circuit and one term of an integrating circuit. Ihe 
coefficient of the integrating circuit turned out to be soro 
indicating this type circuit need not be used. 

In order to show the nature of the feedback function, 5 (p), 
set 3*£(p) = 0 end solve for F Q (p) in equation (30). 

* 0 (p) = (a.i/p+ao+^P+agP 2 ) / L - A-i/P + A Q 4 (41) 



Evaluating L by substituting the values of constants used 
previously, approximating by a power series in (i °° ) and 

neglecting higher order terms gives: 

It s= - 10.59 (i f - 2.33 (i^ ) (42) 



equation (41) can now be evaluated by using the known values of all 
terns, resulting in: 



o 4 

? 0 (p) *s t- 1.06(i~f - 5.84(1^) - s.sol / [ + 10.59(i<->) 

2 

- - 1.13/ (i<0) *«• 2.57 - 12(i^) + 54.77(1^) 



2 . 33 ( 1 ^ )3 
(43) 



therefore, to our order of approximation: 
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A -X = - 1.12 

.Py, = * 2.57 

(44) 

= — 12 

Ag s= + 54.77 

It 13 now necessary to choose a relationship for ? 0 (p) that will 
approach soro for largo co and still retain the desired stabilising 
characteristics for snail u) . Thus the following egression, which 
is arbitrarily chosen, describes a circuit that will stabilise the 
eor.buation ehanber: 

j? 0 (p) = (-1.12/p) (1 + B^p / 1 + Bgp" + Bgp 3 ) (45) 

To evaluate B^, Bg, 3^ in tarns of the constants in equation (44), 
equate equation (45) to (43) giving: 

32^ *= - 2.29 

B « - 10.71 (46) 

b + 23.9 

Thus, 

F 0 (p) « (-1.12/p) (1 - 2.29p / 1 - 10.71p 2 + 23. 9p 3 ) (47) 

This equation for the closed loop circuit gives the desired 
stabilisation to the rocket notor for low frequency pressure 
variations in the burning process which could have caused commotion 
instability. 
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